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ABSTRACT 

We use the most recent Type-la Supernova data in order to study the dark energy - dark matter 
unification approach in the context of the Generalized Chaplygin Gas (GCG) model. Rather 
surprisingly, we find that data allow models with a > 1. We have studied how the GCG 
adjusts flat and non-flat models, and our results show that GCG is consistent with flat case up 
to 68% confidence level. Actually this holds even if one relaxes the flat prior assumption. We 
have also analyzed what one should expect from a future experiment such as SNAP. We find 
that there is a degeneracy between the GCG model and a XCDM model with a phantom-like 
dark energy component. 

Key words: Cosmology :Cosmological Parameters-Observations-Distance Scale-Supernovae 
type la-Method: Data Analysis. 



1 INTRODUCTION 

Recent cosmological observations reveal that the Universe is dom- 
inated by two invisible components. Type-la Supernova observa- 
tions (Riess et al. 1998; Garnavich et al. 1998; Perlmutter et al. 
1999), nucleosynthesis constraints (Buries et al. 2001), Cosmic Mi- 
crowave Background Radiation (CMBR) power spectrum (Balbi et 
al. 2000; de Bernardis et al. 2000; Jaffe et al. 2001), large scale 
structure (Peacock et al. 2001) and, determinations of the matter 
density (Bachall & Fan 1998; Carlberg et al. 1998; Turner 2000) 
allow for a model where the clumpy component that traces mat- 
ter, dark matter, amounts for about 23% of the cosmic energy bud- 
get, while an overall smoothly distributed component, dark energy, 
amounts for approximately 73% of the cosmic energy budget. 

The most interesting feature of this dark energy component 
is that it has a negative pressure and drives the current acceler- 
ated expansion of the Universe (Riess et al. 1998; Garnavich et al. 
1998; Perlmutter et al. 1999). From the theoretical side, great effort 
has been devoted to model dark energy. The most obvious candi- 
date is the vacuum energy, an uncanceled cosmological constant 
[see eg. Bento & Bertolami (1999), Bento et al. (2001)] for which 
lu x = p x / p x — — 1. Another possibility is a dynamical vacuum 
(Bronstein 1933; Bertolami 1986a; Bertolami 1986b; Ozer & Taha 
1987) or quintessence. Quintessence models most often involve a 
single scalar field (Ratra & Peebles 1988a; Ratra & Peebles 1988b; 
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Wetterich 1988; Caldwell et al. 1998; Ferreira & Joyce 1998; Zlatev 
et al. 1999; Binetruy 1999; Kim 1999; Uzan 1999; Amendola 1999; 
Albrecht & Skordis 2000; Bertolami & Martins 2000; Banerjee & 
Pavon 2001a; Banerjee & Pavon 2001b; Sen & Sen 2001 ; Sen et al. 
2001) or two coupled fields (Fujii 2000; Masiero et al. 2000; Bento 
et al. 2002a). In these models, the cosmic coincidence problem, 
that is, why did the dark energy start to dominate the cosmological 
evolution only fairly recently, has no satisfactory solution and some 
fine tuning is required. 

More recently, it has been proposed that the evidence for a 
dark energy component might be explained by a change in the 
equation of state of the background fluid, with an exotic equation 
of state, the generalized Chaplygin gas (GCG) model, rather than 
by a cosmological constant or the dynamics of a scalar field rolling 
down a potential (Kamenshchik et al. 2001; Bilic et al. 2002; Bento 
et al. 2002b). In this proposal, one considers the evolution of the 
equation of state of the background fluid instead of a quintessence 
potential. The striking feature of this model is that it allows for an 
unification of dark energy and dark matter (Bento et al. 2002b). 
Moreover, it is shown that the GCG model may be accommodated 
within the standard structure formation scenario (Bilic et al. 2002; 
Bento et al. 2002b). Concerns about this point have been raised by 
Sandvik et al. (2002), however in this analysis, the effect of baryons 
has not been taken into account, which was shown to be impor- 
tant and allowing compatibility with the 2DF mass power spectrum 
(Beca et al. 2003). Also, the Sandvik et al. (2002) claim was based 
on the linear treatment of perturbations close to the present time, 
thus neglecting any non-linear effects. 

Thus, given it potentialities, the GCG model has been the sub- 
ject of great interest, and various attempts have been made to con- 
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strain its parameters using the available observational data. Studies 
include Supernova data and power spectrum (Avelino et al. 2002), 
age of the Universe and strong lensing statistics (Dev et al. 2002), 
age of the Universe and Supernova data (Makler et al. 2002; Al- 
caniz et al. 2002). The tightest constraints were obtained by Bento 
et al. (2003a) using the CMBR power spectrum measurements from 
BOOMERANG (de Bernardis et al. 2002) and Archeops (Benoit 
et al. 2002), together with the SNe la constraints. It is shown that 
0.74 < A a < 0.85, and a < 0.6, ruling out the pure Chaply- 
gin gas model. From the bound arising from the age of the APM 
08279+5255 source, which is A s > 0.81 (Alcaniz et al. 2002), 
one can get tight constraints, namely 0.81 < A 3 < 0.85, and 
0.2 < a < 0.6, which also rules out the ACDM model. These re- 
sults were in agreement with the WMAP data (Bento et al. 2003b). 
It was also shown that the gravitational lensing statistics from fu- 
ture large surveys together with SN la data from SNAP will be able 
to place interesting constraints the parameters of GCG model (Silva 
& Bertolami 2003). As we shall see in Sections 3 and 4, all these 
constraints are consistent with Supernova data at 95% confidence 
level. 

Recently Choudhury & Padmanabhan (2003) have analyzed 
the supernova data with currently available 194 data points [see 
also Padmanabhan & Choudhury (2003)] and shown that it yields 
relevant constraints on some cosmological parameters. In particu- 
lar, it shows that when one considers the full supernova data set, it 
rules out the decelerating model with significant confidence level. 
They have also shown that one can measure the current value of 
the dark energy equation of state with higher accuracy and the data 
prefers the phantom kind of equation of state, lux < — 1 (Caldwell 
2002). Moreover, the most significant observation of their analy- 
sis is that, without a flat prior, the latest Supernova data also rules 
out the preferred flat ACDM model which is consistent with other 
cosmological observations. In a previous paper, Alam et al. (2003) 
have reconstructed the equation of state of the dark energy compo- 
nent using the same set of Supernova data and found that the dark 
energy evolves rapidly from lo x ~ in the past to a strongly neg- 
ative equation of state (u) x < —1) in the present, suggesting that 
ACDM may not be a good choice for dark energy. More recently, 
other groups have also analyzed these recent Supernova data in the 
context of different cosmological models for dark energy (Gong & 
Duan 2004; Nesseris & Perivolaropoulos 2004). 

In this paper, we analyze the GCG model in the light of the 
latest supernova data (Tonry et al. 2003; Barris et al. 2003). We 
consider both flat and non-flat models. Our analysis shows that the 
problem with the flat model, which has been discussed in Choud- 
hury & Padmanabhan (2003), can be solved in the GCG model in 
a sense that flat GCG model is consistent with the latest Supernova 
data even without a flat prior. We have also analyzed the confidence 
contours for a GCG model, that one expects from a future experi- 
ment such as SNAP. We find that there is a degeneracy between the 
GCG model and a XCDM model with a phantom-like dark energy 
component. 

This paper is organized as follows. In Section 2 we discuss 
various aspects of the generalized Chaplygin gas model and its the- 
oretical underlying assumptions. In Section 3 we describe our best 
fit analysis of the most recent supernova data in the context of gen- 
eralized Chaplygin gas model. Section 4 contains our analysis for 
expected SNAP results. Finally, in Section 5 we present our con- 
clusions. 



2 GENERALIZED CHAPLYGIN GAS MODEL 

The generalized Chaplygin gas (GCG) is characterized by the equa- 
tion of state 
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where A and a are positive constant. For a = 1 the equation of 
state is reduced to so-called Chaplygin gas scenario first studied in 
cosmological context by Kamenshchik et al. (2001). Inserting the 
above equation of state in the energy conservation equation, one 
can integrate it to obtain (Bento et al. 2002b) 
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where p cM) is the present energy density of GCG and ^4,, = 

A/Pir'- 

One of the most striking features of this expression is that, 
the energy density of this GCG, p c h, interpolates between a dust 
dominated phase, p c h oc a~ 3 , in the past and a de-Sitter phase, 
p c h — —pch, at late times. This property makes the GCG model an 
interesting candidate for the unification of dark matter and dark en- 
ergy. Indeed, it can be shown that the GCG model admits inhomo- 
geneities and that, under the Zeldovich approximation, they evolve 
in a qualitatively similar fashion like the ACDM model (Bento et 
al. 2002b). Furthermore, this evolution is controlled by the homo- 
geneous parameters of the model, namely, a and A. 

There are several important aspects of the above equation 
which one should discuss before constraining the relevant parame- 
ters using Supernova data. Firstly, one can see from the above equa- 
tion that A s must lie in the range < A s < 1. For A s = 0, GCG 
behaves always as matter whereas for A s = 1, it behaves always 
as a cosmological constant. Hence to use it as a unified candidate 
for dark matter and dark energy one has to exclude these two pos- 
sibilities resulting the range for A s as < A s < 1. 

To have an idea about the possible range for a, one has to 
consider the propagation of sound through this fluid. Given any 
Lagrangian C(X, <f>) for a field <f>, where X — ^g 11 " '4>,it4>, v , the 
effective speed of sound entering the equations for the evolution of 
small fluctuations is given by 

C,x 

p,x C,x + 2XC,xx 
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Thus, for a standard scalar field model which has canonical kinetic 
energy term like £ = X — V(<f>), the speed of sound is always equal 
to 1 irrespective of the equation of state. But for a Lagrangian con- 
taining a non-canonical kinetic energy term, one can have a sound 
speed quite different from 1. Actually, even > 1 is possible 
which physically means that the perturbations of the background 
fluid can travel faster than light as measure in the preferred frame 
where the background is homogeneous. For a time dependent back- 
ground field, it does not lead to any violation of causality as the 
underlying theory is manifestly Lorentz invariant (Erickson et al. 
2002). 

For GCG it has been shown that the equation of state (1) can 
be obtained from a generalized version of the Born-Infeld action 
(Bento et al. 2002b) 



C = -A 



lKl+a) [l-W v <t>,A, v t +a)/2a ] a/{1+a) , 



(4) 



which for a = 1 leads to the Born-Infeld action. If one computes 
c 2 3 for this action one can get using Eq. (2) the present value as 
Cso = a A s . As As is always positive, its restricts a to only positive 
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Figure 1. Confidence contours in the a — A s parameter space for flat uni- 
fied GCG model. The solid and dashed lines represent the 68% and 95% 
confidence regions, respectively. The best fit value used for M is -0.033. 



values. In all previous work, a has been restricted to a value up to 
1. But one can see from the above expression for c 2 that as < 
A s < 1, the maximum allowed value for a can be surely greater 
than 1 and that also depends on the value of A s , e.g for A 3 = 
0.5, the allowed range for a is < a < 2. Notice also that the 
dominant energy condition p + p > is always valid in this case. 
Furthermore, there is no big rip in the future and asymptotically the 
Universe goes toward a de-Sitter phase. 

Hence, on general grounds, restricting a up to 1 is not a very 
justified assumption. Moreover, this restriction arises mainly by 
considering the present day value of cj which is not an important 
epoch for structure formation. In general c 2 in this model is a time 
dependent quantity and in such cases it is not very proper to con- 
strain a with the present day value of c 2 . 

There is another reason for not restricting a up to 1 . It has been 
shown by Kamenshchik et al. (2001) that one can also model the 
Chaplygin gas with a minimally coupled scalar field with canoni- 
cal kinetic energy term in the Lagrangian density. Performing this 
exercise for the GCG, leads to a potential for this scalar field of the 
form 

V = Vc ) e 3(Q - 1) ^cosh(^) 2 /^+ 1 ) +C osh(^r 2 ^ +1 >] (5) 

where Vb is a constant and m — 3(a + 1). For a = 1, one recovers 
the potential obtained by Kamenshchik et al. (2001). Now as we 
have discussed earlier, for a minimally coupled scalar field with a 
canonical kinetic energy term, the value of c 2 . is always 1 irrespec- 
tive of the equation of state. Hence if one considers this kind of 
scalar field to model GCG, there is no such restriction on a coming 
from the sound speed. 

In what follows, we shall consider that A s lies in the range 
< A s < 1 and the only constrain on a that we shall consider is 
that it takes positive values. We should also point out that the a — 
case corresponds to the ACDM model. 

The Friedmann equation for a non-flat unified GCG model in 
general is given by 

H 2 = Hi [Q ch (As + (1 - + zf (1+a) ) 

+ Q k (l + z) 2 ] , (6) 
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Figure 2. Same as Figure 1, but with a wider range for a. 



where Ho is the present day value of the Hubble constant. Q c h 
and fit are the present day density parameters for GCG and the 
curvature. For a flat Universe Q k = and Q ch = 1, whereas for 
the non-flat case Q. c h = 1 — fifc. 



3 RECENT SUPERNOVA DATA AND THE BEST FIT 
ANALYSIS 

To perform the best fit analysis of our GCG model with the recent 
Supernova data, we follow the method discussed by Choudhury & 
Padmanabhan (2003). As far as the Supernova observation is con- 
cerned, the cosmologically relevant quantity is the apparent mag- 
nitude, m, given by 

m(z)=M + 5\og w D L (z) , (7) 

where Dl = —<1l{z), is the dimensionless luminosity dis- 
tance, and the luminosity distance cIl(z) is given by cIl(z) = 
r(z)(l + z) where r(z) is the comoving distance. Also M = 

M + 51og 10 ( xmIc ) + 25 wnere M is the absolute magnitude 
for the Supernova which is believed to be constant for all Type-la 
Supernova. 

In our analysis, we take the 230 data points listed in Tonry et 
al. (2003) along with the more recent 23 points from Barris et al. 
(2003). Also, as discussed by Choudhury & Padmanabhan (2003), 
for low redshifts, data might be affected by the peculiar motions, 
making the measurements of the cosmological redshifts uncertain; 
hence we shall consider only those points with redshifts z > 0.01. 
Moreover, since it is difficult to be sure about the host galaxy ex- 
tinction, A v , we do not consider points which have A v > 0.5. 
Hence in our final analysis, we consider only 194 points, which are 
similar to those considered by Choudhury & Padmanabhan (2003). 

The Supernova data points given by Tonry et al. (2003) and 
Barris et al. (2003) are listed in terms of luminosity distance 
log 10 dh{z) together with the corresponding error a\ ogl0 d L (z). 
These distances are obtained assuming some value of A4 which 
may not be the true value. Hence, in our analysis we shall keep it 
as free parameter while fitting the data. 

The best fit model is obtained by minimizing the quantity 
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Using best fit value A = 0.89 
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Figure 3. Confidence contours in the — a parameter space for non flat 
unified GCG model. The solid and dashed lines represent the 68% and 95% 
confidence regions, respectively. The best fit value used for M is -0.033. 
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where M = M — Mobs is a free parameter denoting the differ- 
ence between the actual M and the assumed value Mobs in the 
data. As discussed by Choudhury & Padmanabhan (2003), we have 
also taken into account the uncertainty arising because of the pecu- 
liar motion at low redshift by adding an uncertainty of Av = 500 
kins" 1 to a? oglo dL (z), 



0k>g lc d L (z) -» cr? ogw dL (z) + 



1 Aw 
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(9) 



This correction is more effective at low redshifts, i.e for small val- 
ues of dL- 

In our subsequent best fit analysis, the minimization of (8) is 
done with respect to M ,a, A 3 and fife. The parameter M is a 
model independent parameter and hence its best fit value should 
not depend on a specific model. We have checked that when min- 
imizing (8) with respect to M , the best fit value for M for all 
of the models considered here is —0.033 which is also consistent 
with that obtained by Choudhury & Padmanabhan (2003). Hence, 
in our subsequent analysis, we shall use always this best fit value, 
M = -0.033. 



3.1 Flat case 

For this case, we assume Qk = and consider only two parame- 
ters, a and A s . We first restrict a to be < 1. In Figure 1, we have 
shown the 68% and 95% confidence contours in a — A s parameter 
space. The best fit values for [a, A s ] are given by [0.999, 0.79]. 
The best fit value of a is very close to its upper limit since the 
actual best fit value lies in the region beyond a = 1. Also, up to 
68% confidence level, the a = 0, i.e. the ACDM case, is excluded 
although it is consistent at 95% confidence level. 

Next we allow a to vary for a wider range. The Figure 2 repre- 
sents the same as Figure 1 but with a taking a wider range. The best 
fit values for [a, A 3 ] are now [3.75, 0.936]. This high value of a 
might just be a statistical fluke though, as the confidence regions ex- 



Figure 4. Same as Figure 3, but in the Q k — A s parameter space. 



hibit a very shallow valley along the a direction. Here, also at 68% 
confidence level, the a — 0, ACDM case, is excluded, although it 
is consistent at 95% confidence level. It should be noted that both 
Choudhury & Padmanabhan (2003) and Tonry et al. (2003) also 
found some conflict between the ACDM model and the SN la data, 
namely, that when imposing a flat Universe prior, the data ruled out 
the vacuum energy as an allowed dark energy component, and ac- 
tually favoured a phantom energy component. Here we see that the 
GCG model fits the data well, and, as mentioned previously, with- 
out the theoretical complications that plague the phantom energy 
model, namely the dominant energy condition is not broken, and 
there is no big rip singularity in the future. 

It is also clear from the minimum value for \ 2 obtained in 
these two cases, that when one allows to vary a beyond 1, one 
obtains a better fit to the Supernova data. 



3.2 Non-flat case 

Another problem that the ACDM has with the new SN la data is 
that without a flat prior, a flat ACDM Universe has also been ruled 
out at 68% confidence (Choudhury & Padmanabhan 2003). Tonry 
et al. (2003) argued that this was probably due to some overlooked 
systematic error, since a small systematic error of 0.04 mag was 
able make the flat ACDM consistent with the data. Here we attempt 
to find if the GCG model might alleviate this problem. 

For this, we allow a non-vanishing curvature in our model. 
We now have three parameters in our model namely, a, A 3 and the 
density parameter for the curvature at present, fife. First we assume 
that our Universe deviates slightly from the flat model assuming Qk 
to vary between [-0.1,0.1]. In this case the best fit values for a, A s 
and n fe are [2.87, 0.89, -0.099]. It suggests that the data prefers a 
negative curvature. In Figure 3, we have shown the 68% and 95% 
confidence contours in the Qk-a plane assuming the best fit value 
for A s , whereas in Figure 4, we have shown the same contours in 
the £lk-A s plane assuming the best fit value for a. 

Both a and A s are constrained significantly and 68% confi- 
dence limits on a and A s are [1.6 3.625)] and [0.856 0.946)] re- 
spectively. It shows that the data prefers higher value for a and the 
ACDM limit (a = 0) is excluded, both for 68% and 95% confi- 
dence limit. Also one can see from both figures that the flat case 
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Figure 5. Confidence contours in the fij, — a parameter space for non-flat unified GCG model. The solid and dashed lines represent the 68% and 95% 
confidence regions, respectively. The best fit value used for M is -0.033. Figure (a) is for the best fit mean value of A s , whereas (b) and (c) are for the values 
in the wings of 68% confidence limit. 




Figure 6. Same as Figure 5, but in the Qfc — A s parameter space. Figure (a) is for the best fit value of a whereas (b) and (c) are for the values in the wings of 
68% confidence limit. 



fife = is consistent with the data for both 68% and 95% confi- 
dence level, but for a higher value of a. 

We allow now more curvature in our model and consider the 
range for to vary as [-1, 1]. In this case the best fit values for 
a, A s and Slk are [0.73, 0.65, -0.999] and the resulting Xmin is 
197.99 showing that there is an improvement in the quality of fit. 
The 68% confidence limits on a and A s are [0.052 1.056] and [0.62 
0.81], respectively. Also the allowed range for a shifts more sig- 
nificantly towards smaller values and data do allow the model to 
get closer to ACDM (a = 0.052) as well as to the Chaplygin gas 
model (a = 1). 

In Figure 5, we have shown the 68% and 95% confidence con- 
tours in fife -a plane assuming the best fit value for A s , as well as 
for its values in the wings of 68% confidence limit. It shows that 
the allowed range is quite sensitive to the parameter A s . For smaller 
values of A s (but within its 68% confidence limit) the flat model 
(fife = 0) is more inconsistent with the data and negative curvature 
is preferred. But for higher values of A s , e.g A s = 0.81 which 



still falls within the 68% confidence limit, it shows that flat mod- 
els as well as models with small but both positive and negative 
curvature, are allowed [consistent with WMAP bound on Sltotai, 
0.96 < Cltotal < 1-08 (Spergel et al. 2003)] but for non zero val- 
ues of a. It suggests that assuming a GCG model, one can alleviate 
the problem of consistency with a flat universe as pointed out by 
Choudhury & Padmanabhan (2003) for a ACDM model. 

In Figure 6, we have shown the same contours but now in the 
Q,k-A 3 plane assuming the best values for a as well as its values at 
the wings of the 68% confidence limit. Figure (6c) is for a = 0.052 
which is almost a ACDM model, and it is quite similar to what ob- 
tained by Choudhury & Padmanabhan (2003) for a ACDM model. 
It shows that model that behaves more like ACDM (a = 0.052) 
is inconsistent with a flat universe at 68% confidence level. On the 
other hand, model that deviates more from the ACDM model, Fig- 
ure (6a) and (6b), is consistent with the flat universe with better 
confidence level. Like Figure 5, it again shows that even if one does 
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Table 1. SNAP specifications for a two year period of observations. 



Redshift Interval 


Number of SNe 


0.0-0.2 


50 


0.2-1.2 


1800 


1.2-1.4 


50 


1.4-1.7 


15 
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Figure 7. Expected confidence regions for SNAP for a GCG fiducial model 
with 1 — A s = 0.25 and a = 1. The solid and dashed lines represents 
the 68% and 95% confidence regions respectively. The left contours are for 
GCG and the right ones are for XCDM. For GCG, the parameter space is 
1 — A 3 Vs. a whereas for XCDM, it is Sl m Vs. fi. We have marginalized 
over M. 



not take a flat prior, unlike the ACDM model, flat GCG model is 
consistent with the supernova data up to 68% confidence level. 



4 EXPECTED SNAP CONFIDENCE REGIONS 
4.1 Method 

To find the expected precision of a future experiment such as SNAP, 
one must assume a fiducial model, and then simulate the experi- 
ment assuming it as a reference model.This allows for estimates of 
the precision that the experiment might reach [see Silva & Berto- 
lami (2003) and references therein for a more detailed description 
of the method employed here]. Let us then assume a fiducial model 
and functions \ 2 based on hypothetical magnitude measurements 
at the various redshifts. In this case, 

1 2 



X 2 (model) = ^ 



[ m model( 2 ')- m fid( z ' 
a\z) 



(10) 



where the sum is made over all redshift bins and m(z) is as speci- 
fied in Section 3. 

Here together with the GCG model we also analyze a flat 
XCDM model (that is a Cold dark matter model together with a 
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Figure 8. Same as figure 7 but for a XCDM fiducial model with Q m 
0.49 and (3 = -oj x = 1.55. 



dark energy with an equation of state p x — u> x p x ). We aim to show 
that if the Universe is indeed described by the GCG model, the 
fitting of a XCDM model to the data will reveal that the cosmic 
expansion is drawn by a phantom-like dark energy component. 

To fully determine the \ 2 functions, the error estimates for 
SNAP must be defined. Following Weller & Albrecht (2002), we 
assume that the systematic errors for the apparent magnitude, m, 
are given by 

°- 02 

Osys = yg-Z , (11) 

which are measured in magnitudes such that at z = 1.5 the sys- 
tematic error is 0.02 mag, while the statistical errors for m are es- 
timated to be Osta = 0.15 mag. We place the supernovae in bins of 
width Az « 0.05. We add both kinds of errors quadratically 



,(%) = 



^sys {%i 



+ 



" sta 



(12) 



where m is the number of supernovae in the i'th redshift bin. The 
distribution of supernovae in each redshift bin is, as before, taken 
from Weller & Albrecht (2002), and shown in Table I. 

Summarizing, for each fiducial model, the method used, con- 
sists in the following 

(i) Choose a fiducial model. 

(ii) Fit the XCDM model to the mock data, and obtain the re- 
spective confidence regions. 

(iii) Repeat the previous step to the GCG. 



4.2 Results 

In Figures 7 to 9 we show the confidence contours for the GCG 
and XCDM model for future SNAP observation taking different 
fiducial models. In all these Figures j3 = — us x . Also, along x-axis 
in all of these figures, we have plotted 1 — A s (instead of A s as 
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Figure 9. Same as figure 7, but for a ACDM fiducial model with f2 m 
0.30. 



1 — A s represents Q. m when a = 0) for the GCG and f2 m for the 
XCDM. 

As mentioned above, our main aim is to explicitly show that 
a GCG Universe might appear as a XCDM Universe with a dark 
energy component that has a phantom-type equation of state. To do 
so, we have considered two fiducial models. The first corresponds 
to a Chaplygin model (a = 1), with 1 — As — 0.25. If one attempts 
to fit a XCDM model to the data (Figure 7), one finds that data 
favour a larger amount of matter than expected and a phantom-type 
dark energy component. This is fully consistent with Figure 13 of 
Tonry et al. (2003). 

To further examine the degeneracy between models, in Fig- 
ure 8 we have repeated the procedure assuming a XCDM fiducial 
model, with Q. m = 0.49 and (3 = 1.55 (w = —1.55). From ex- 
amining Figures 7 and 8, one can see that both models appear es- 
sentially identical. Also for GCG, the confidence regions for both 
fiducial models are quite identical to what we have shown earlier in 
Figure 2 for the current Supernova data. 

In Figure 9 we used a fiducial ACDM model, with Q m = 0.3. 
As can be seen, the confidence regions are completely different 
from those found in the two previous cases. Also, for the GCG 
model the confidence regions are quite different from what we have 
earlier in Figure 2 for the current Supernova data, further hinting 
that indeed the ACDM model is not a good description of the Uni- 
verse. 

To illustrate the degeneracy between the GCG model and the 
XCDM model with a phantom like equation of state (lj x < —1), 
we Taylor expand the luminosity distance as, 



z + 



(1 



go - 3qo + jo)z 3 + ... 



(13) 



where go is the deceleration parameter related with the second 
derivative of the scale-factor and jo is the so-callled jerk parameter 
(Visser 2003) related with the third derivative of the scale-factor. 
This is also one of the statefinder variables r proposed in Sahni et 
al. (2002). The subscript "0" means that quantities are evaluated at 



present. The jerk parameter is related with the deceleration param- 
eter qo as 



„ 2 dq I 
jo = q + 2q + — |o 



(14) 



For the GCG model, one can calculate qo and |o to get 



(1-A S )-1 
A 3 (l-A a )(l + a), 



GCG 

Qo 

dq , GCG 

whereas for XCDM model they turn out to be 

§< 

9 



(15) 



DE 



-(l - u x (n m - 1)) - 1 



dq.DE 
d~z l ° 



(16) 



For the previous Supernova data obtained by Perlmutter et al. 
(1999) and Riess et al. (1998) for low redshifts (z < 1), it is suf- 
ficient to consider the first two terms in the series expansion of the 
luminosity distance cLl{z) given above. In that case, one can see 
from the expression of qo for GCG, that Supernova data can only 
constrain A s , as go, for the GCG model, is independent of a. More- 
over, in order to have degeneracy between the GCG and XCDM 
model, the go parameter of these two models must be equal, which 
results that: 



A s = u> x (0, m — 1) 



(17) 



Thus, if one gets a bound on A s by fitting GCG model with low 
redshift Supernova data, then the same data can be fitted by a host 
of XCDM models (including ACDM model) provided the above 
equation is satisfied. Hence, for low redshifts (z < 1), GCG model 
is degenerate with all kind of different dark energy models with 
constant equation of state, including the ACDM model. 

Now, as one goes to higher redshifts, which is the case for the 
current data that we are considering in this paper, one also has to 
consider the higher order terms in the series expansion of cLl (z) . As 
far the data we are studying in our paper, it is enough to consider 
terms up to order z 3 in the series expansion of cLl. Hence, in order 
to have a degeneracy between GCG and XCDM even for the high 
redshifts, the jerk parameter jo also has to be equal for the two 
models, which effectively means -p|o has to be equal for the two 
models. Using this together with the equation (17), one finds that 



-a(l- A s ) - 1 
(l + a)(l- A s ) 



(18) 



l + a(l-A s ) 

In the above equation, lo x and Q, m are the equation of state 
and the density parameter of the XCDM model which is degenerate 
with a GCG model with the corresponding a and A s parameters, 
for higher redshifts. Now one can see that for any GCG model (a > 
0), the corresponding equation state has to be always phantom type 
{ijj x < —1) as < A s < 1. This shows that although for the low 
redshift data, GCG model is degenerate with all kinds of constant 
equation of state dark energy model including the ACDM model, 
for higher redshift, the GCG is degenerate with a XCDM model but 
only with a phantom type of equation of state. 

Finally, aiming to further illustrate the degeneracy of the GCG 
model with the XCDM model for large red-shifts, we plot in Fig- 
ure 10 the behaviour of the dimensionless luminosity distance Dl 
as function of redshift for four different best fit models, namely: 
ACDM (w x = -1, fi m = 0.3), XCDM (w x = -2.07 Q m = 
0.51), Chaplygin model (a = 1, A s = 0.77) and GCG (a = 3.75, 




0.5 1 1.5 



z 

Figure 10. Dimensionless luminosity distance Dl /z as function of redshift 
for four different best fit models. The solid line is for ACDM. The dotted, 
dashed and dashed-dot lines are for Chaplygin, XCDM and GCG model, 
respectively. The values for the parameters of these models used in the plot 
are given in the main text. 



A s — 0.936) respectively. We have actually plotted the streched 
Dl/z as function of redshift in order to graphically show the de- 
generacy. 



5 CONCLUSIONS 

We have analyzed the currently available 194 supernova data points 
within the framework of the generalized Chaplygin gas model, re- 
garding GCG as a unified candidate for dark matter and dark en- 
ergy. We have considered both, flat and non-flat cases, and used the 
best fit value for M = —0.033 which is independent of a specific 
model, throughout our analysis. For the first time, we have crossed 
the a = 1 limit for the GCG model and try to see whether the data 
actually allow it or not. 

For the flat case, we have studied both cases, restricting a to 
the range < a < 1 and also without any restriction on the upper 
limit q. From Figures 1 and 2, it is quite clear that data favours 
a > 1, although there is a strong degeneracy in a. Also the quality 
of fit improves substantially as one relaxes the a = 1 restriction. 
Moreover, the minimum values allowed for a and A s at 68% con- 
fidence level are [0.78, 0.778], which excludes the a = 0, ACDM 
case, although there is no constrain on a at 95% confidence level. 

Moreover, if one does not assume a flat prior for the analy- 
sis, our study shows that the flat GCG model for a values suffi- 
ciently different from zero, is consistent with the Supernova data 
up to 68% confidence level. It also allows small, both positive and 
negative curvature, making the GCG a somewhat better description 
than the ACDM model. This is consistent with recent result which 
shows that without a flat prior, a flat ACDM model, which is oth- 
erwise consistent with different cosmological observations, is not a 
good fit to the supernova data (Choudhury & Padmanabhan 2003). 
Moreover, the fact that GCG is a better fit to the Supernova data 
than ACDM, is consistent with the result of Alam et al. (2003), 



who also have reconstructed a similar kind of evolving equation of 
state for the dark energy from the latest Supernova data. 

We have also studied the confidence contours for a GCG and 
XCDM model expected from the future SNAP observation assum- 
ing different fiducial universes. In this regard, the degeneracy be- 
tween the GCG model and a phantom-like dark energy scenario has 
made obvious in Section 4, where we have shown that when fitting 
a XCDM model to a GCG universe, the data will favour a phan- 
tom energy component, and vice- versa. This degeneracy is also il- 
lustrated analytically through the expression for the luminosity dis- 
tance cLl as function of redshift. We have shown that for higher red- 
shifts, GCG model is completely degenerate with a XCDM model 
with a phantom type of constant equation of state. We mention that 
it has already been noted in Maor et al. (2002) that time varying 
equations of state might be confused with phantom energy, and 
here we show that this is true for the GCG, without breaking the 
dominant energy condition and without a big rip singularity in the 
future. It should also be noted that with the exception of a cosmo- 
logical constant, most dark energy models predict a time varying 
equation of state, therefore a constant dark energy equation of state 
might not be the best parametrization for dark energy. We have also 
shown that for a ACDM fiducial model, confidence regions for a 
GCG model, which are expected from future SNAP experiment, 
are quite different from what we have shown in Section 3.1, sug- 
gesting that SN la data does not favour the ACDM model. 

Thus, our study shows that the generalized Chaplygin gas 
model is a very good fit to the latest Supernova data both with or 
without a flat prior. With future data, one expects the error bars to be 
reduced considerably, but we still expect that Supernova data will 
favour a generalized Chaplygin gas model with high confidence. 
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